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Abstract 

Analytical expressions for width and conductance peak distributions for quan- 
tum dots with multi-channel leads in the Coulomb blockade regime are pre- 
sented for both limits of conserved and broken time-reversal symmetry. The 
results are valid for any number of non-equivalent and correlated channels, 
and the distributions are expressed in terms of the channel correlation matrix 
M in each lead. The matrix M is also given in closed form. A chaotic billiard 
is used as a model to test numerically the theoretical predictions. 
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Advances in nanostructures technology make possible the manufacture of semiconduc- 
tor devices known as quantum dots J!]] where electrons are confined to very small two- 
dimensional regions. By connecting external leads to such devices it is possible to study 
their electronic transport properties: the conductance can be measured as a function of the 
Fermi energy and/or as a function of a magnetic flux through the dot. Of particular interest 
are dots that are weakly coupled to leads, due to the presence of barriers at the interface 
between dot and leads. In these cases the resonance widths T are small compared to their 
mean spacing A. At low temperatures kT < A, only one quasi-bound level participates in 
the conduction process. The resulting Coulomb blockade peaks of the conductance || are 
equally spaced but their amplitude exhibits order of magnitude fluctuations. 

For dot sizes that are smaller than the electron-impurity mean free path, conductance 
fluctuations are determined by the dot geometry. We discuss dot geometries which display 
classical chaotic motion. In such cases one can model their transport properties using the 
concepts of chaotic scattering ||. Our results should also hold for weakly disordered dots 
in the quasi-zero-dimension limit. In Ref. || a statistical description of the conductance 
peaks in the Coulomb blockade regime was developed using the i?-matrix formalism |5j 
and assuming that the dot's wavefunctions are described by random matrix theory (RMT). 
Conductance and decay width distributions were derived for one-channel leads, as well as for 
symmetric leads with several equivalent and uncorrelated channels. However, a more general 
theory would encompass leads with an arbitrary number of correlated and non-equivalent 
channels. Progress in this direction was made in Ref. [[?], where by modelling leads as 
point contacts || spatial wave function correlations were taken into account. Using the 
supersymmetry method ||, exact formulas for the width and conductance peak distributions 
were obtained for two-point contact leads. However, the derivation was restricted to the case 
of broken time reversal symmetry, and could not be applied to leads with finite width or 
with more than two-point contacts. On the experimental side, we note that conductance 
distributions are becoming accessible; for ballistic open dots (r ^> A) such distributions 
were recently measured || , and similar experiments are underway in the Coloumb blockade 
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regime. 

In this letter we present exact formulas for the width and conductance peak distributions 
for leads with any number of channels that are in general correlated and non-equivalent. 
These distributions are obtained both for conserved and broken time-reversal symmetry and 
are completely characterized in terms of the channel correlation matrix M. Our results are 
valid for both the point-like contacts and the continuous extended leads models. We are 
able to treat the most general case because our methods are based exclusively on RMT, 
which is technically simpler than the supersymmetry method used in Ref. 0. To test our 
theory we use a chaotic billiard, the Africa ||10|| , for which the statistical distributions of 



one-channel leads were recently studied in detail (TT|. This model is particularly useful to 
study systems with strong channel correlations. In contrast, the correlations between nearest 
points in the discretized Anderson model of a disordered dot were too weak to produce 
any significant change in the width distribution of two-point leads (as compared with the 
uncorrelated channels distribution). We remark that since the partial width is analogous to 
the wavefunction intensity (see Eq. (^) below), our results for the partial and total width 
distributions can be directly tested in the microwave cavity experiments |T2|, where the 



wavefunctions intensities are measured at several points and are spatially correlated. 

Provided that T <C kT < A, which is typical of experiments fl|||, the conductance peak 
amplitude for a two lead geometry is given by 

G = — —— g , where g = , A A , (1) 
h 2kT y ' y + ' v ; 

and (r^) is the partial decay width of a resonance A into channels of the left (right) lead. 
Each lead can support open channels so that T 1 ^ = J2 C |7ca^| 2 > where 7^ is the partial 
amplitude to decay into channel c on the left (right). i?-matrix theory gives || 

7cA=(^) V 7^ $ :( r )^(r), (2) 

where ^\ is the resonance wave function inside the dot (scattering region), $ c is the wave 
function of an open channel c in the lead (asymptotic region), and the integral is taken over 



the contact boundary between the lead and the dot. P c and k c are the channel penetration 
factor to tunnel through the barrier and the longitudinal wave number, respectively. An 
alternative modelling assumes that the quantum dot is connected to the leads by one or 
more point-like contacts ||. Every such point contact r c is considered as a channel and the 
corresponding partial amplitude is |7|] 

7c A = (« c ^A/vr) 1/2 vI/ A (r c ), (3) 

where A is the area of the dot, A is the mean resonance spacing and a c is the coupling 
parameter of the point contact to the dot . Expanding a resonance wave function with energy 
e in a fixed basis of states with that energy ^\(r) = J2\ VVP/x( r ) (the sum is truncated to N 
terms, typically much larger than A), the partial width to decay to channel c can be expressed 
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as a scalar product 7 cA = (</> c |V>a) = E^^VVj where 0*^ = (h 2 k c P c /rnj JdS $*(r)p M (r) 

1 /2 

in the i?-matrix formalism and 0* M = (a c AA/n) p^(r c ) in the point-like contact model. 
Expressing the width as a scalar product allows us to treat the extended lead and the 
point-like contact models in an equivalent manner. 

The resonance states ^\ are assumed to have GOE- or GUE-like statistical properties, 
depending on the symmetry class to which the dynamics in the dot corresponds [Q. This 
assumption is valid both for dots with chaotic dynamics and for weakly disordered dots. 
The eigenvector components (ipi,ip2, ■ ■ ■ > iPn) = V* (i n the following we shall omit the eigen- 
vector label A) are therefore randomly distributed P(iJ)) oc 5(J2^=i IVvJ 2 — -0 EH- ^ ne j° m t 
distribution of the partial width amplitudes 7 = (71, 72, ... , 7a) for A channels is given by 

p( 7 ) oc /dm 5 (jr i^/vi 2 - 1 ) n % c - (0 C |V>» , (4) 

\n=l J c=l 

where the metric is D[if)] = H^dip^ for the GOE and D[if)) = Y\^ =1 dip*dip fl /2'Ki for the 
GUE. To evaluate (ffl) we transform cf) c = Y,c> 4>c'Fc'c to obtain a new set of orthonormal 



channels (0 C |0 C ,) = 5 cc i. In the limit iV — > 00, provided that A<iV, one finds [|1J 

P( 7 ) = (detM)- /3 /V^ tM -^, (5) 
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where M = (NF^F) 1 . The distribution ([5]) is normalized with the measure D[~y] = 
nli ^7c/2tt for the GOE (f3 = 1) and £>[ 7 ] = Uc=i dj*d^ c /2m for the GUE ((3 = 2). Note 
that for both ensembles the joint partial width amplitudes distribution is Gaussian, and it 
follows that M is just the correlation matrix of the partial widths M cc / = 7*7c' = jj(<p c \<p c r). 
In general the channels are correlated (non-orthogonal) and non-equivalent, i.e. have differ- 
ent average partial widths. 

Recalling Eqs. (§) and @, the spatial autocorrelation function C(Ar) = 



\[ / *(r)\I/(r + Ar) /|\I/(r) | 2 plays a central role in deriving explicit expressions for the cor- 
relation matrix M. For fully chaotic systems with time reversal symmetry, Berry obtained 
C(Ar) semiclassically, assuming that classical orbits cover uniformly the energy surface 
||16|| . For an eigenstate of a chaotic billiard with energy s = h 2 k 2 /2m, C(Ar) = Jo(k\ Ar\). 
Alternatively, a resonance at energy e can be expanded inside the dot in the fixed basis 
p M (r) = A" 1 ^ 2 exp(ik^ • r), where the /i's correspond to N different (random) orientations 
of k. In the spirit of RMT, assuming that the expansion coefficients ip^ are Gaussian (this 



hypothesis was confirmed for the Africa billiard |T7[) and using ip*^^ = iV" -1 ^/, one can 
rederive Berry's result. See also Ref. [18]. The presence of a small magnetic flux $ in- 
troduces corrections to C(Ar) which are small in the semiclassical limit and are of order 
h 2 ($ / &o) 2 j (2mAe) -C 1. Thus, the correlation matrix is given by 

M cd = a2(fccP ^ Pc,)1/2 JdsJdS'K(r)Jo(k\ r - r '\)Mr') (6) 
for the finite width leads, and 

M cc , = A(acac ' )1/2 J (k\r c -r' c \) (7) 

7T 

for the point contacts model. 

We turn next to the calculation of the total width distribution P(T) in a given lead 
that supports A channels and is characterized by a correlation matrix M. Using (|5|) and 
T = S c | Tc 1 2 5 the characteristic function of P(T) is readily obtained and we find 

j roo P~^~ 

per) = — dt • (8) 

2tt7-oo [det(J - 2itM/(3)f 2 



The matrix M is hermitian and positive definite, so that its eigenvalues w 2 are all positive. 
Since T is invariant under orthogonal (unitary) transformations, P(T) depends only on w 2 
and Eq. (H) can be evaluated by contour integration. When all eigenvalues of M are non- 
degenerate, we find for the GUE case 



-r/«. 



TT(--- " 

C '}c W d W c 



(9) 



For two equivalent channels (M n = M 22 = T/2), Eq. (|9|) coincides with the result of 



PauEit) = ^ e --/(H/n sinh (^Lf) , (10) 

where T = T/T, and / = M 12 / ^/M^M^ measures the degree of correlation between the two 
channels. 

For the GOE case, the integral in @ is still straightforward but cannot be expressed 
in a simple form as above. Labelling the inverse eigenvalues of M in increasing order 
Wi 2 < < . . ., the contour integral gives 



c / m=l " ( 2w 2m-i) 1 

(11) 



n^fV-i)iT) 
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2w 2) VV s '=2m\2w^ i T ) 



where, for an odd number of channels A, we define l/2u^ +1 — > oo. For the case of two 
equivalent but correlated channels, Eq. (]TT|) reduces to 

Pco Mt) = -^e-^I (ft), (12) 



where / = Mu/^/MuM^ and Jo is the modified Bessel function of order zero. 

To test the RMT predictions, we modelled a quantum dot by the Africa billiard fllQ| , |TT|j . 
The shape of this billiard is determined by the image of the circle of radius one in the complex 
2-plane under the conformal mapping w(z) — (z + bz 2 + ce lS z 3 )/Vl + 2b 2 + 3c 2 . We studied 
the case b = 0.2, c = 0.2 and 5 = n/2 which has a classical fully chaotic phase space [ID]. For 



the statistical study of the case of broken time reversal symmetry, we consider the billiard 



threaded by an Aharonov-Bohm flux line $ = a$o [ 1Q| , |TT | . We choose a = 1/4 where the 



resonance fluctuations are known to be GUE-like |TTJ. The study considers 50 eigenstates 
starting from the 100£h. 

To investigate the eigenfunction amplitude correlations C(Ar) we have used the eigen- 
functions of the billiard with von Neumann boundary conditions. The results are shown in 
the inset to Fig. p] where the correlations in the model (solid line) compare well with the 
theoretical result Jo{kAr) (dashed line). In what follows, we assume for simplicity that the 
penetration factors k c P c are channel- independent. We first studied extended leads by taking 
the contact region of the lead and the dot to have a finite width D. In this case, the channels 
are defined by the allowed transverse momenta ttc/D (c = 1, . . . , A) where A = intffc-D/V]. 
To guarantee that the correlation matrix M is the same for different eigenfunctions of the 
billiard, we choose D such that kD is constant and scale the partial amplitude (0) by k. We 
find that the channels in this case are weakly correlated. Thus, the total width distribution 
is similar to the case of uncorrelated equivalent channels, which gives a \ 2 distribution with 
A (GOE) or 2A (GUE) degrees of freedom. However this changes if the barrier penetration 
factors have a strong energy dependence. Our model calculations agree nicely with the RMT 
predictions (see Fig. [I]). 

We also studied the model of leads with A point-like contacts by choosing for each lead 
a sequence of A equally spaced points on the Africa boundary. Thus, according to (0), M is 
completely determined by A and fc|Ar| (where |Ar| is the distance between two neighboring 
points). In Fig. |2| we compare the results of this model (histograms) with the theoretical 
predictions (solid lines) for A = 4 and different values of fc|Ar|. 

To calculate the conductance distribution, which is the measurable quantity for quantum 
dots, we assume that the left and right leads are uncorrelated and characterized by M l and 
M r , respectively. We then use P(g) = J dT l dV r 5(g - TT^ /(T l + T r )) P(T l )P(T r ), where 



P(r) is given by ([|) or (|TTJ) . In the absence of time reversal symmetry we find 
Pgue(9) = — T2^e «= - 

{YlcVcYldWd) C A 



nd-i)n(i-i 

c'^c 



w- 



2g 

V c W d 



1 / V c W d . , , 

2 \w d v c 



V c W d 



(13) 



where v c (wd) are the eigenvalues of M l (M r ) and Kq (Ki) is the modified Bessel functions 
of order zero (one). 

The result of [§]|| is a special case of ( |13"1) for one channel leads with = T r (i.e. 
V\ = w±), while the distribution of Ref. is obtained for two (equivalent) channels leads. 
For time reversal symmetric systems, we also obtained a closed formula for P(g), which has 
similar structure to (|13D. Fig. [5] shows a comparison between the theoretical conductance 
distributions and those calculated for the Africa billiard for A-point symmetric leads with 
fc|Ar| = 1 and for various values of A. 

In conclusion, we have derived in closed form the width and conductance peak distri- 
butions in a quantum dot, for leads with any number of correlated and/or non-equivalent 
channels, and in the presence or absence of time reversal symmetry. The only required input 
to determine the distributions is the channel correlation matrix M, for which an explicit ex- 
pression was obtained. Our results for the decay widths could also be applied to compound 
nucleus reactions in the limit of isolated resonances, where M is evaluated by the optical 
model. 

We acknowledge A.D. Stone and H. Bruus for useful discussions and for the use of 
their Africa billiard computer program and E. R. Mucciolo for discussions. This work was 
supported in part by DOE Grant DE-FG02-91ER40608 and by NSF. C.H.L. thanks the 
Center for Theoretical Physics at Yale University for its hospitality. 
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FIGURES 

FIG. 1. Total width distributions P(T) for extended leads with A = 4. T is measured in 
units of its average value. The solid lines are the theoretical distributions and the histograms are 
results from the Africa billiard for the (a) GOE and (b) GUE cases (see text). The insert shows 
the spatial wavefunction correlation C(Ar) calculated for the Africa (solid line) compared with 
Jo(k\Ar\) (dashed line). 

FIG. 2. Total width distributions P(F) for A = 4 point-contact leads, (a) GOE, fc|Ar| = 0.5; 
(b) GUE, jfe|Ar| = 0.5; (c) GOE, fc|Ar| = 1 and (d) GUE, fc|Ar| = 1. The solid lines correspond to 
the theoretical distributions, the dashed lines to uncorrelated channels and the histograms to the 
Africa billiard calculations. 

FIG. 3. Conductance peak distributions P(g) for A point-contact symmetric leads (M l = M r ) 
with fc|Ar| = 1. (a) GOE, A = 2; (b) GUE, A = 2; (c) GOE, A = 4 and (d) GUE, A = 4. The 
convention for the lines is as in Fig. 2 . 
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